Nonlinear current helicity fluxes in turbulent dynamos and alpha quenching 
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Large scale dynamos produce small scale current helicity as a waste product that quenches the large scale 
dynamo process (alpha effect). This quenching can be catastrophic (i.e. intensify with magnetic Reynolds 
number) unless one has fluxes of small scale magnetic (or current) helicity out of the system. We derive the 
form of helicity fluxes in turbulent dynamos, taking also into account the nonlinear effects of Lorentz forces due 
to fluctuating fields. We confirm the form of an earlier derived magnetic helicity flux term, and also show that it 
is not renormalized by the small scale magnetic field, just like turbulent diffusion. Additional nonlinear fluxes 
are identified, which are driven by the anisotropic and antisymmetric parts of the magnetic correlations. These 
could provide further ways for turbulent dynamos to transport out small scale magnetic helicity, so as to avoid 
catastrophic quenching. 
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Large scale magnetic fields in astrophysical bodies are 
thought to be generated by dynamo action involving helical 
turbulence and rotational shear A particularly impor- 

tant driver of the mean field dynamo (MFD) is the a effect, 
which in the kinematic regime is proportional to the kinetic 
helicity of the turbulence. A question of considerable debate 
is how the a effect gets modified due to the backreaction of 
the generated mean and fluctuating fields? It is especially 
important to understand whether a suffers catastrophic (i.e. 
R m -dependent) quenching, since R m , the magnetic Reynolds 
number, is expected to be typically very large in astrophysical 
systems. Recent progress has come from realizing the impor- 
tance of magnetic helicity conservation in constraining this 
nonlinear saturation 01 01- 

Recall that in the MFD theory, one splits the magnetic field 
B into a mean magnetic field B and a small scale field b = 
B B, and derives the mean-field dynamo equation [ 1 ] 

fyo 

— = vx (UxB + e-nJ). (i) 



Here £ = u x b is the turbulent electromotive force (emf), 
J = V x B/(j,q the mean current density, /io the vacuum 
permeability (assumed unity throughout the rest of the paper), 
r/ the microscopic resistivity and the velocity U = U + u 
has also been split into mean U and small scale turbulent u = 
U—U velocities. Finding an expression for the correlator £ in 
terms of the mean fields is a standard closure problem which is 
at the heart of mean field theory. In the two-scale approach 1 2] 
one assumes that £ can be expanded in powers of the gradients 
of the mean magnetic field. For isotropic helical turbulence, 
this gives £ k, aB — rjtJ, where in the kinematic limit, a — 
«k = —\tu) ■ it, proportional to the kinetic helicity (u> = 
V x u) and rjt — \i~u 2 is the turbulent magnetic diffusivity 
proportional to the specific kinetic energy of the turbulence, 
with r being the velocity correlation time. 

The kinematic theory has to be modified to take account 



of the backreaction due to the Lorentz forces associated with 
the generated large and small scale fields. Treatments of the 
back -reaction have typically used the quasi-linear approxima- 
tion or closure schemes to derive corrections to the mean-field 
dynamo coefficients (cf. |4j,|5ll and |2[ for a review). It is then 
found that the a effect gets "renormalized" by the addition of 
a term proportional to the current helicity, j ■ b, of the small 
scale fields; that is a = qk + au, where om = {T~/3po)j ■ b 
|3l • (Here po is the density of the fluid and j the small scale 
current density.) At the same time there is no modification to 
rjt at lowest order |4]. 

In order to constrain otyi, in simulations in a periodic box 
or in systems that involve no boundaries, it has proved useful 
to take recourse to the evolution equation for the small scale 
magnetic helicity h — (a ■ b), where a is the vector potential 
and ( } denotes volume averaging. We have in such situations, 
dh/dt = -2(5-5) - 27/ (j • b). In the stationary limit, this 
predicts (£ ■ B) = —r){j ■ b), which tends to zero as r\ — > 0, 
for any reasonable spectrum of current helicity. This leads to a 
catastrophic quenching of the turbulent emf parallel to B. Of 
course while evolving to this stationary state, some B will be 
generated. But its value, at the end of the kinematic regime, 
still turns out to be small, if the scale separation is large |6]. 
It has been suggested that such quenching can be avoided if 
the system has open boundaries and is inhomogeneous, since 
one could then have a flux of small scale helicity out of the 
system which helps maintain a non-zero (£ ■ B) |2l lia uS- 
It is therefore important to calculate such fluxes in a general 
manner, taking into account also the effect of Lorentz forces. 
This is the aim of the present work. 

Indeed, Vishniac and Cho 01 derived an interesting flux 
of helicity, which arises even for nonhelical but anisotropic 
turbulence. We derive a generalized form of the Vishniac- 
Cho flux (henceforth VC flux) in the evolution equation for 
j bio include also nonlinear effects of the Lorentz force and 
helicity in the fluid turbulence. As we see, the VC flux can 



2 



also be thought of as a generalized anisotropic turbulent dif- 
fusion. Further, due to nonlinear effects, other helicity flux 
contributions arise generated by the anisotropic and antisym- 
metric part of the magnetic correlations. 

One immediate problem with previous approaches was that 
in open systems with boundaries, h is not gauge invariant and 
one has to consider instead the gauge-invariant relative mag- 
netic helicity, say Iir, defined by subtracting the helicity of 
a reference vacuum field 1 1 0] . The flux of relative helicity 
is cumbersome to work with for arbitrarily shaped bound- 
aries. Also the concept of a density of relative helicity is not 
meaningful, since fiR is defined only as a volume integral. 
In order to avoid these problems it is advantageous to con- 
sider instead the current helicity j ■ b and its flux. The cur- 
rent helicity density and its flux are directly gauge invariant, 
locally well defined, and are in fact the observationally mea- 
sured quantities in say the solar context. Furthermore it isj b 
which directly enters the expression for the nonlinear a effect 
Also for isotropic small scale fields the spectra of small 
scale current (Cfe) and magnetic helicities (Hk) are related by 
Hk = Ck/k 2 . For these reasons we consider here directly the 
current helicity evolution and use the current helicity flux as a 
'proxy' for the magnetic helicity flux. 

Current helicity evolution. Consider the evolution of the small 
scale current helicity, j ■ b = eijkbidjbk, which is explic- 
itly gauge invariant. We assume that the correlation tensor 
of fluctuating quantities (u and b) vary slowly on the sys- 
tem scale, say R. Consider the equal time, ensemble average 
of the product f(xi)g(x2). The common dependence of / 
and g on t is assumed and will not explicitly be stated. Let 
f(ki) and 5(^2) be the Fourier transforms of / and g, re- 
spectively. We can express this correlation as f(xi)g(x2) = 
J $(/, g, k, R) e ik r d 3 k, with 

$(f,g,k,R) = J f(k + \K)g(-k + \K) e lKR d 3 K. 

(2) 

Here we have defined the difference r = X\ — x 2 and the 
mean R — \{x\ + x^), keeping in mind that all two-point 
correlations will vary rapidly with r but slowly with R 1 1 1 ] . 
Also k = \{k\ — ^2) and K = k\ + k2- In what follows we 
require the correlation tensors, Vij (k,R) = iij,k,R), 
of the u field, rriij (k, R) = &(bi, bj, k, R) of the b field and 
the cross correlation Xjk(k,R) = &(uj,bk,k, R) between 
these two fields, in Fourier space. The turbulent emf is given 
by £i{R) = e ijk J Xjk(k,R)d 3 k. 

In order to compute the current helicity evolution of 
dj(x) -b(x)/dt, we use the induction equation for b in 
Fourier space, 



db k {k) 
dt 



= Ckpqtqlmikp J Uj (fe - k')B m (k') d fc 

+G k (k)~ v k 2 b k {k). 



(3) 



Here, G = Vx(ux6-«xb)is the nonlinear term. (We 
also neglect the velocity shear due to U compared to that due 



to it.) A tedious but straightforward calculation gives J5I 

d f r _ _ 

j b = eijk %Xlk kj(k ■ B) — XifcVj(ife • B) 



df 



-ikj B ■ V\ik + 2ikjXpkV P Bi 



d 3 fc + T c .(4) 



We have written out explicitly only that part of the helicity 
evolution driven by the coupling of the turbulent emf to the 
mean magnetic field. This is because we are particularly inter- 
ested in turbulent helicity fluxes driven by an inhomogeneous 
£ and B. The Tc term represents the triple correlations of the 
small scale u and b fields and the microscopic diffusion terms 
that one gets on using Eq. yjl- The handling of the triple cor- 
relations needs a closure approximation. But we will not need 
to explicitly evaluate these terms to identify the helicity fluxes 
we are interested in; i.e. those which couple £ and B, and so 
continue to write this term as Tc- 

In order to calculate the current helicity evolution, using 
Eq. @, we have to calculate also xik- This has been done 
in detail by id. (One adopts a closure approximation 
whereby triple correlations, TJ^, which arise now in the evolu- 
tion equation for dxik/dt, are assumed to provide relaxation 
of the turbulent emf or xik and one takes Tik — ~Xik/T, 
where t is the relaxation time (cf. also 1 13])). We concentrate 
below on nonrotating but helical turbulence. For such turbu- 
lence we have from 1 3 , El , Xik = r\ik, where \ik is given by; 
see Eq. (10.30) in f^T 

kk = —tic ■ B(vik -mik) + \B - V(vik + mtk) 



+Bi s m s 



Bk, s Vls 



h R 



m, ^ra.s 



dvi k | dmi k \ 
dk s J 



Ok 



2 ^2 s ,pTTlpk • 



(5) 



We use this in what follows. 

Let us denote the four terms under the integral in Eq. @ by 
Ai, A-2, A3 and A4, respectively. In A\, due to the presence of 
eijk, only the antisymmetric parts of the tensors vik and mih 
survive, and these are denoted by vfj, and mf^., respectively. 
Also note that the last term above vanishes because it involves 
the product eijkhkj = 0. All the other terms of Eq. @ al- 
ready have one R derivative, and so one only needs to retain 
the term in xik = 7"l;fe which does not contain R derivatives. 
In ^3 one can then use the fact that V ■ B = to write it as a 
total divergence. We now turn to specific cases. 
Isotropic, helical, nonrotating turbulence. Let us first recon- 
sider the simple case of isotropic, helical, nonrotating, and 
weakly inhomogeneous turbulence. For such turbulence, the 
form of the velocity and magnetic correlation tensors is given 
by lHH [l2il . In evaluating the fc-integrals, only terms which 
involve integration over an even number of ki survive. Also, 
in terms which already involve one Ri derivative, one needs 
to keep only the homogeneous terms in and mik- Fur- 
ther in the presence of eijk all terms symmetric in any pair of 
the indices vanish. Taking account of these considerations, it 
turns out that only the homogeneous part of the velocity and 
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magnetic correlation tensors given in survive. The 

homogeneous part of these correlations is 



gives 



k^ kj 



E(k,R)- € -^F(k,R), (6) 



and a similar expression for my with functions say M(k, R) 
and N(k, R) replacing E and F, respectively. Here Ank 2 E 
and 47rfc 2 M are the kinetic and magnetic energy spectra, re- 
spectively, and Airk 2 F and 4irk 2 N are the corresponding he- 
licity spectra. They obey the relations u 2 = 2 J Ed 3 k, 
u-V xu = 2f Fd 3 k, b 2 = 2 J Md 3 k, and j~b = 
2 J Nd 3 k. With these simplifications we have, after carry- 
ing out the angular integrals over the unit vectors fc, = fcj/fe, 



+B p B mtj / k m k p {v^ k +m^)d 3 fc 

d 3 fc. (10) 



/ MfcWfc + 



All the other terms, A2, A3 and A4, cannot be further simpli- 
fied. They are explicitly given by 

A 2 = -Tei 3k B p B s j / k,kp{vfc - mt k ) d 3 k, (11) 



M = ±B 2 



k 2 (F~N) d 3 k + %B ■ J 



k 2 {M + E) d 6 k. 

(7) 

In the case of isotropic turbulence, the second and third 
terms, A2 and A3, are zero because, to leading order in R 
derivatives, the integrands determining A2 and A3 have an 
odd number (3) of fc, 's. The fourth term is given by 



A 4 = § J • Bt I k 2 [E - M] d 3 k. 



(8) 



Adding all the contributions, A\ + A2 + A3 + A4, we get for 
the isotropic, helical, weakly inhomogeneous turbulence, 



\B t I k 2 {F-N)d 3 k+\j-BT I k 2 Ed 3 k+T c . 



d_- 

at' 

(9) 

We see that there is a nonlinear correction due to the small 
scale helical part of the magnetic correlation to the term oc 
B . But the nonlinear correction to the term oc J ■ B has 
canceled out, just as there is no such correction to turbulent 
diffusion 0], 

As pointed out above, in the isotropic case the magnetic he- 
licity spectrum is Hk = Gk/k 2 . So the first two terms of the 
current helicity evolution equation Eq. l|9} can be interpreted 
as representing the effects of exactly the source term —2£ B 
which is obtained for the magnetic helicity evolution. Also for 
this isotropic, but weakly inhomogeneous case, one sees that 
there is no flux which explicitly depends on the mean mag- 
netic field. 

Anisotropic turbulence. Let us now consider anisotropic tur- 
bulence. In the first term A\ in Eq. l|4}, one cannot now as- 
sume the isotropic form for the velocity and magnetic correla- 
tions. But again, due to the presence of ey^, only the antisym- 
metric parts of the tensors vik and mik survive. Also the last 
term in Eq. does not contribute to A\ because it involves 
the product eijkfakj = 0. One can further simplify the term 
involving k derivatives by integrating it by parts. Straightfor- 
ward algebra, and a judicious combination of the terms then 



As 



Tti jk B p B s / kjk p (v lk - m lk ) 



d 3 fc, (12) 



A4 = 2Te ljk B p Bi iS / kjkp(v sk - m sk )d 3 k. 



(13) 



Adding all the contributions, A\ + A2 + A3 + A4, we get 

d 



01 



j b = 2e jlk T 



B P B S I ikjk p k s (vik - m^) d 3 k 



2B p B k , s I k J kpvf s d 3 k - B p B s j I kgkpm^d^k 



V s I B P B S I kjkpmt k ) d 6 k 



(14) 



Here vf s — \{vi s + v s i) is the symmetric part of the velocity 
correlation function. 

Let us discuss the various effects contained in Eq. H41 for 
current helicity evolution. The first term in Eq. dl4> repre- 
sents the anisotropic version of helicity generation due to the 
full nonlinear a effect. In fact, for isotropic turbulence it ex- 
actly will match the first term in Eq. l|9}. The second term in 
Eq. J14i gives the effects on helicity evolution due to a gener- 
alized anisotropic turbulent diffusion. This is the term which 
contains the VC flux. To see this, rewrite this term as 



dj ■ b 



Of 



ATe 3 ikBpB k , s I kjk p vf s d 3 k 



V T v +4TB k e k ijB p , s I k 3 k p vf s d 3 k. (15) 



Here the first term is the VC flux, T s = <p S pkB p B k , where 
4> sp k is a new turbulent transport tensor with 



i> sp k = -^rejik J kjkpvf s d 3 k = -ATUj k \7 p u s 



(16) 



Obviously, only components of <j> spk symmetric in p and k en- 
ters in the flux T s . The second term in Eq. d!5l > is the effect 
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on helicity due to 'anisotropic turbulent diffusion'. (We have 
not included the large scale derivative of vi s to the leading or- 
der.) Strictly speaking, T v is a current helicity flux, but if 
we define the spectrum of the magnetic helicity flux by divid- 
ing the spectrum of the current helicity flux by a k 2 factor, 
Eq. d!6l > for T v leads exactly to the magnetic helicity flux 
given in Eqs. (18) and (20) of Vishniac and Cho |7|. 

This split into helicity flux and anisotropic diffusion may 
seem arbitrary; some support for its usefulness comes from 
the fact that, for isotropic turbulence, T v vanishes, while the 
second term exactly matches with the corresponding helicity 
generation due to turbulent diffusion, i.e. the J ■ B term in 
Eq. Of course, we could have just retained the non-split 
expression in Eq. dl5b . which can then be looked at as an effect 
of anisotropic turbulent diffusion on helicity evolution. Also, 
interestingly, there is no nonlinear correction to the VC flux 
from the small scale magnetic field in the form of, say, a term 
proportional to mf s ; just as previously, there was no nonlinear 
correction to turbulent diffusion in lowest order! 

Finally, Eq. H41 also contains terms (the last two) involv- 
ing only the antisymmetric parts of the magnetic correlations. 
These terms vanish for isotropic turbulence, but contribute 
to helicity evolution for nonisotropic turbulence. The last 
term gives a purely magnetic contribution to the helicity flux, 
but one that depends only on the antisymmetric part of mik- 
Note that such magnetic correlations, even if initially small, 
may spontaneously develop due to the kinematic a effect or 
anisotropic turbulent diffusion and may again provide a he- 
licity flux. More work is needed to understand this last flux 
term better fl4l[l5ll . Preliminary simulations of helical turbu- 
lence with shear and open boundaries suggest that the sign of 
the VC flux agrees with that of the small scale current helicity 
flux, but that its magnitude may only account for about 25% 
of the actual flux 1 16]. The existence of the VC flux has also 
been verified in simulations of nonhelically driven shear flow 
turbulence 1 17], but its magnitude was too small to produce 
dynamo action. 

In conclusion, we have derived helicity fluxes in turbulent 
dynamos, taking also into account the nonlinear effects of 
Lorentz forces due to the fluctuating field. To avoid gauge 
ambiguities, we have followed the current helicity evolution. 
We confirm the form of the helicity flux found by Vishniac and 
Cho, who used the first order smoothing approximation. We 
note however that it is more correctly interpreted as a current 
helicity flux and not as a flux of relative magnetic helicity. In 
addition we have found that the corresponding turbulent co- 
efficient does not get renormalized due to nonlinear effects, 
just as is the case of turbulent diffusion. Additional nonlinear 
fluxes have been identified as being driven by the anisotropic 
and antisymmetric parts of the magnetic correlations. These 
could provide further ways for turbulent dynamos to transport 
out small scale magnetic helicity so as to avoid catastrophic 



i? m -dependent quenching. It remains to calculate these fluxes 
in specific circumstances and also verify their presence in di- 
rect numerical simulations of turbulent dynamos. 
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